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Abstract 

We study the partition function of the free Sp{N) conformal field theory recently 
conjectured to be dual to asymptotically de Sitter higher-spin gravity in four-dimensions. 
We compute the partition function of this CFT on a round sphere as a function of a 
finite mass deformation, on a squashed sphere as a function of the squashing parame- 
ter, and on an S'^ x geometry as a function of the relative size of S*^ and S^. We 
V/~) find that the partition function is divergent at large negative mass in the first case, and 

^ for small 5*^ in the third case. It is globally peaked at zero squashing in the second 

O case. Through the duality this partition function contains information about the wave 

^ function of the universe. We show that the divergence at small occurs also in Ein- 

^ stein gravity if certain complex solutions are included, but the divergence in the mass 

^ parameter is new. We suggest an interpretation for this divergence as indicating an 

5^ instability of de Sitter space in higher spin gravity, consistent with general arguments 

that de Sitter space cannot be stable in quantum gravity. 
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1 Introduction 

The dS/CFT correspondence is perhaps the most well-known program for incorporating 
cosmologically relevant spacetimes into some form of holographic theory [1, 2, 3, 4, 5].^ One 
proposal for the dS/CFT dictionary-^ [5] is that the Hartle-Hawking wave function [17] in the 
bulk is equal to the partition function of a Euclidean CFT deformed by various operators 
with finite coefficients (see also [18, 19, 20, 21, 22, 23, 24, 25]). In this paper we study that 
proposal in a recent explicit realization of the dS4/CFT3 correspondence [26], in which the 
bulk is Vasiliev's higher spin gravity [27, 28] in dS4 [29, 30]. We explicitly compute the 
partition function with some finite deformations turned on and discuss the implications for 
the wave function of Vasiliev's universe. 

The duality of [26] was inspired by similar efforts in the context of AdS4 higher spin 
gravity [31, 32, 33, 34, 35, 36, 37]. The dual CFT is a (Euclidean) free theory of anti- 
commuting scalar fields x"? with fiat space Lagrangian [38, 39, 40] 

CcFT = Inatda^d'x' . (1.1) 
iPor other proposals see [6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. 

^For a recent review on some aspects of asymptotically de Sitter spaces see [16]. 
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Here Qat is the symplectic form, with a, 6 G {1, 2, ... , A^}, and the x"" transform as an Sp{N) 
vector.'^ Repeated indices are summed over. N is related to the bulk cosmological constant 
Atuik = ^/ils as AT ~ l/{GAbuik)- 

To get a bulk spectrum that matches Vasiliev's, the operator content is restricted to the 
Sp{N) singlet sector. This is built from a set of even-spin bilinear currents of the general 
form QabX"'d^-^ where the derivatives can act either to the left or the right and 

are symmetrized with various traces removed such that for s > 2 the current is symmetric, 
traceless, and conserved. We will be especially interested in the spin zero 'current' 

= -^nabx^x' , (1.2) 

which has conformal dimension one, and the spin two current Tij, which is the conserved 
traceless stress tensor [41]^ 

Tij = ^Qab [-3d,x''d,x' + X'^d^d.x' + ^^Ax^d^'] . (1.3) 

This stress tensor has conformal dimension three. 

J^^'^ is dual to a bulk scalar field with mass rn^uik^'ds ~ 2- Tij is dual to the bulk 
graviton. More generally, the spin s current is dual to a massless spin s field in the bulk. In 
relating the bulk and boundary there is a subtlety for the scalar in that the dimension of J° 
is less than 3/2, so what is usually called the "alternate quantization" in AdS/CFT is in play 
here. Turning on a field theory source a for J° does not correspond to Dirichlet boundary 
conditions for the bulk scalar field at the future boundary; as we explain in appendix A 
it instead corresponds to fixing a linear combination of and its canonical momentum 11 at 
the late-time cutoff surface. 

To be explicit, recall the dS/CFT dictionary in [5]: 

where ^ hh ['Pi 9] is a Wheeler-de Witt wave function for the scalar field configuration at the 
boundary of spacetime to be and the induced metric on the boundary to be g computed 
from an integral over compact bulk geometries with no other boundaries. The equality really 
is meant asymptotically as the "cutoff" parameter e goes to zero. A is the dimension of the 
■^Note that in our conventions N must be even. 

^These operators are renormalized in a way we will make explicit below. 
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dual operator O, and Soft is understood to include divergent local terms in the sources 
a and g, which are tuned to match the divergent local pieces of "^hh- ^ In the alternate 
quantization instead of fixing the value of (p at the boundary we instead fix a quantity which 
for the Sp{N) model is 

- e^n = ^TiVeV . (1.5) 

Here 11 = — is the momentum conjugate to (j) acting on the wave function. In Fefferman- 
Graham gauge, where the metric approaches ds"^ = {ids/TY {—dT^ + Qijdx'^dx^) as T — > 0^, 
we have 

n = T-29t0 . (1.6) 

The coefficients of (f) and 11 in (1.5) are chosen to ensure conformal invariance with the 
correct dimension for the dual operator O. We can thus interpret the free Sp{N) partition 
function as the same Hartle-Hawking wave function projected onto a basis of eigenstates of 
the Hermitian operator — e'^U. Unlike in AdS/CFT there is only one bulk quantization; 
changing the scalar boundary conditions is simply changing the quantum mechanical basis 
for the wave function. For more details on this point see appendix A. 

In AdS/CFT it is usually sufficient to restrict attention to the partition function with 
arbitrary infinitesimal sources turned on for the boundary operators; for us this would have 
the form 




where . . . are linear couplings to the higher spin currents and also any counterterms. 

This linearized expression allows computation of bulk correlators with fixed future bound- 
ary conditions a = 0, hij = [19, 24, 25], but in order to make sense of the full wave function 
proposal 

^HH [^V, e^'^Qij, ...] = ZcftW, 9ij, • • •] , (1-8) 

we need to understand how to generalize (1.7) to finite values of a, hij, etc. In this paper we 
will consider only nonzero cr and hij , since these couple to relevant and marginal operators in 
the Wilsonian sense. Intuitively a gives a position-dependent mass term and hij introduces 
a background metric in the CFT. The higher spin currents are irrelevant, and turning them 

^It is somewhat unsatisfying that the scheme-dependent local terms in the field theory partition function 
must be determined by comparison to the bulk, but the fact is that the local terms in the wave function 
do contain physical information. From the field theory point of view there is not any particular reason to 
choose one scheme over another, with the exception of trying to subtract the divergences altogether. This 
choice of scheme is incorrect here, since the divergences really do exist on the bulk side. One might hope 
that it is possible to formulate a purely field theoretic condition which determines the "correct" coefficients 
for divergent local terms in any particular cutoff scheme. 
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on at finite strength introduces considerable UV subtlety into the definition of the partition 
function. The presence of higher-spin symmetry may constrain these operators in a strong 
enough way that the partition function remains well-defined despite their irrelevance, but 
we will not attempt to demonstrate this here. For the most part we will just set all of the 
higher-spin sources to zero.^ 

Our interpretation of this wave function will be the same as that of Hartle and Hawking 
[17], which is that its square computes conditional probabilities. So for example 



gives the probability that the induced metric on some spatial surface is gij and that the 
scalar and its canonical momentum obey (1-5), given that the higher spin fields are all zero. 

The outline for the remainder of the paper is as follows: in section 2 we introduce a 
renormalization procedure for dealing with finite sources and use it to compute the partition 
function on a round as a function of a constant source a for J^^\ We discover that it 
diverges at large negative a. In section 3 we compute the partition function of the Sp{N) 
model on a squashed three sphere as a function of the squashing. We find that it is globally 
peaked at zero squashing. In section 4 we compute the partition function on S'^ x S"^ as a 
function of the relative size of S"^ and S^, finding that it diverges at small S^. In section 
5 we recall the bulk wave function of a free massive scalar and observe that there is no 
divergence of the type found in section 2. We then do an Einstein gravity dS version of 
the Hawking-Page calculation of the wave function on S'^ x S^, finding that we are able to 
produce a divergence of the same type as we found from the dual field theory in section 4. In 
section 6 we return to the scalar divergence; we show that it implies that the critical version 
of the Sp{N) model at finite N does not exist on S^. We interpret this as an instability of 
higher spin gravity, modulo a certain technical point about integrating over the higher-spin 
fields that we are unable to decisively resolve. In section 7 we briefiy point out a relationship 
between the Sp{N) model on quotients of "Hs and bulk 'big bang' cosmologies. Section 8 
is our conclusion, and in appendix A we give a detailed treatment of double-trace fiow in 
dS/CFT and the dS interpretation of the "alternate quantization". 

point that has been emphasized by Steve Shenker is that single-trace irrelevant operators do not 
appear in conjectural duals to more conventional dS theories that do not have higher spin massless fields in 
the bulk. In particular the operator dual to a scalar with positive mass squared never has a real part which 
is greater than d, and the operators dual to gauge fields and the graviton are always at most marginal. 
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2 Finite Deformations and the Wave Function on the Three Sphere 



The obvious guess for how to turn on finite hij and a is to just put the theory (1.1) on a 
background Qij and turn on a position dependent mass a. Including the conformal mass 
term this gives 



Share 9ij ) X] 



ab 



(2.1) 



Justifying this expression however involves some care. In particular putting the theory on 
a background metric gij involves nonlinear terms in the perturbation hij\ for example the 
conformal mass is required to preserve tracelessness of the stress tensor. One might worry 
that nonlinear terms in a similarly need to be included in the bare action of the field theory 
in order to properly reproduce the bulk physics, or even that terms involving higher powers 
of the curvature tensor and its derivatives need to be included.' Fortunately in the free 
Sp{N) model any term involving o"^ is always irrelevant in the Wilsonian sense, as are most 
things involving the curvature tensor. The only relevant (or marginal) terms that can be 
added depend only on the sources: 

Sloe = AA^ j d^x^ j d^^Va^la] + ^^J ^^^v^^' (2-2) 

where A is the energy cutoff. We are here assuming that the theory is cutoff in a geometric 
way so that all terms are diffeomorphism invariant. The linear term in a has a simple 
interpretation as a freedom to shift the identity part of the composite operator x^, and the 
other terms are renormalizations of the energy momentum tensor. There is also a relevant 
operator (J^^^y and a marginal operator {J^^^y, but these must have zero coefficient to stay 
in the free theory in the IR when a = 0. 

The need to impose the singlet constraint also leads to an additional subtlety, which is 
typically dealt with by coupling the theory to Chern-Simons theory and taking the limit k — )■ 
oo. Fortunately this does not affect the partition function on simply-connected manifolds 
like and S^, since there are no nontrivial flat connections, so we will postpone further 
discussion of the Chern-Simons sector until section 4.*^ 



^We thank Dan Jafferis and Douglas Stanford for discussions of this issue. 

®A point that has been emphasized to us by Simeon HeUerman is that on the pure Chern-Simons 
partition function vanishes in the hmit k — > oo. Since we should really multiply (2.3) below by this partition 
function, the result becomes trivial. For this reason we have to imagine studying the theory at large but 
finite fc, as emphasized in [42]. This vanishing however cancels out of any relative probabilities computed 
on the topological three-sphere. In this paper we will always compare different metrics and scalar field 
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In the cases where we may neglect the Chern-Simons sector, we may formally evaluate 
the path integral: 

ZcFT[<J,gij\ = e-^'"^ J I5;^e-^''— = e-^'-det (^ ~^ ttV^t^^^^^ ) • (2-3) 
Here we have defined the path integral measure as 

n 

where c° are the Grassman coefficients of n-th eigenfunction of the conformal Laplacian 
— + ^R[g] in the mode expansion of x", and a is a dimensionless coefficient which can be 
adjusted by rescaling x aiid is thus a possible wavefunction renormalization. 

The coefficients of the local terms in Sioc are in principle determined by matching to a 
bulk calculation of the wave function, but in Vasiliev gravity this is beyond reach since the 
action is not known. We will instead observe that for typical bulk theories the divergences 
in '$HH[4>,g] are pure phase, so that they cancel in the probability \E'\E'*. For example we'll 
see this is true for pure Einstein gravity or a free bulk scalar in section 5 below. The general 
reason is that if the bulk Lorentzian action is real, the dominant contribution to the action 
at large spatial boundary comes from real classical field configurations near the boundary 
which produce a phase e*'^. Without the action we cannot be completely sure that this 
works in Vasiliev theory, but in any case it is required for there to be a well-defined de Sitter 
invariant probability distribution I^E^/^hP at late times. From here on we will assume that the 
divergent local pieces in '^hh[(P, g] are indeed pure phase. We can thus define new partition 
function Z finite, which is related to Zcft by a divergent local phase, by 

/ — + cr + -R\q] \ 

ZfiniteW, gij] = e-^=* det ( ' j . (2.5) 

Here is of the same form as (2.2), but with A, B, C, and a now chosen to cancel the 
UV divergences of the determinant. This new definition is extremely convenient since it can 



configurations on the same topology, either the three sphere or 5^ x 5^, so our results for these probabilities 
will be excellent approximations to the "true" finite but large k theory. 

^This can also be justified by analytic continuation from the AdS version of Vasiliev's theory if we assume 
that that all local terms are real in a Euclidean AdS calculation of the partition function as a function of 
boundary sources. By dimensional analysis the coefficients of these terms are always an odd power of £ads, 
divided by G, and multiplied by a power series in G£~^^. Continuing iads i^ds then gives local terms 
which are purely imaginary for real sources. 
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be unambiguously implemented on the field theory side, and with our phase assumption we 
have 

l^/f/fl^ = \ZcFt\'^ = \Zfinite\^ (2.6) 

SO we may still extract the probability distribution from Z finite- 
2.1 Partition Function with a Constant Mass 

We first discuss the simple case of a constant mass deformation cto on a round S'^. We will 
take a simple hard cutoff in the sum over modes: 

21ogZ,,„.,, = -2Set + AT l)Mog (^^^^4^^^) • (2-7) 

We have used that the spectrum of the Laplacian on 5"^ is + 2) for £ = 0, 1, 2, . . ., that 
the degeneracy is (£ + 1)^, and have chosen units where the sphere radius is 1. To ensure that 
the cutoff is geometric, we insist that it is really on eigenvalues of the conformal Laplacian: 

3 ^ 

^maxi^max ~^ 2) ~l~ ^ ^ ■ (2-8) 

First consider the case where ctq = 0. The sum can be done in Mathematica, and its 
asymptotic expansion for large A is 

2N~hogZ finite = - 2N'^Set - ^(1 + 31oga)A=^ - log a - ^(1 + 71oga)A 

+ ^ (^log4 - 21oga - + 0(1/A). (2.9) 



Cancellation of the quadratic divergence requires a = 1, and cancellation of the cubic and 

187r2 ^^^^ -° - 1887r2 



linear divergences fixes A = —tP-2 and B = — 2 - We thus find 



N A . 3C(3) \ 



log Zf,nUe = ^ [log 4 - , (2.10) 

which is the continuation — )■ — of the well-known result in the 0{N) model [43, 44]. If 
we had not defined A in a geometric way as in (2.8) we would have gotten the wrong finite 
part. 
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Figure 1: A plot of the probability distribution \ Zfinite\^ = \^hh\'^ as a function of a constant mass 
perturbation a on the 5^, for N = 2. 

Turning the mass ctq back on and imposing the cutoff as 

3 

Cax(Cax + 2) + - + do = A^ (2.11) 

a similar argument sets (7 = ^ and gives 

2,ogZ,„,. ^ f (log 4 -my^p,,,r [log (l + j^j^) - 



O"0 



+ 2) + ,, 



3 

(2.12) 

The sum on the right-hand side can be done by differentiating with respect to ctq, performing 
the sum, and then integrating. The constant can be fixed by matching to (2.10) when ao = 0. 
The final result is 



1 7 ^ 

log Z Ur. 



^finite 



487r2 



67r2(l - 4(To) log - e-^^Vi-4<^oj + 12Li3 (^e^^^^T^^ 
+ zvrVl - 4^0 (vr'(l - 4(To) + 12Li2 {e-'-^/^^^y^ ] . (2.13) 



Here Li2 and Lis are polylogarithms. We plot the square of this result in figure 1. We will 
make a few comments, and return to the bulk interpretation later: 

• There is a local maximum at a = 0, which is the conformally invariant point. Zero- 
mode fluctuations are suppressed, which suggests that the dS solution of higher-spin 
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gravity is perturbatively stable. In fact since the two-point function of J^^^ is negative- 
definite this perturbative stabihty extends to the other modes of a as welL Using a 
concavity argument analogous to one we will outline in our discussion of the squashed 
sphere below, we can even prove that this maximum is a global maximum over all 
functions cr(x) which are everywhere greater than — |. 

• The wave function squared vanishes for a = + 2) — |; this is because at these 
values there are fermion zero modes so the determinant vanishes. 

• At large constant positive a we have Z ~ e-^f'^^ so the probability distribution is 
exponentially suppressed. We won't show it here, but this behavior is just that of the 
partition function on M^; this makes sense since when the mass is large the fluctuations 
should not be sensitive to the radius of curvature. 

• At large constant negative a the wave function is divergent; the peaks between the zeros 
grow exponentially in a like \Z\ ~ e"^''^°s2. This suggests that the wave function is 
non-normalizeable in the negative a direction, although we will not be able to decisively 
say this without a discussion of integration in other field directions that we elaborate 
on below. We nonetheless see that the peak of the probability distribution is not at 
the dS invariant point. We are inclined to interpret this as an instability of dS space 
in higher-spin gravity; the equilibrium field distribution at late times is far from its 
dS value. We discuss this point further in section 6. It is worth mentioning that the 
growth does not start until 0(1) values of a, which our boundary conditions (1-5) tell 
us are actually bulk field values that are 0{\/N). Thus from a bulk point of view the 
instability is a finite effect. For comparison we will show in section 5 below that 
there is no such divergence in the wave function for a free bulk scalar in de Sitter space. 

3 The squashed partition function of the Sp{N) model 

In this section we will compute the partition function of the Sp{N) model on a squashed 5"^, 
with (7 = 0. This is a purely metrical deformation, the space remains topologically an 
and we can continue to neglect the Chern-Simons sector. 

^'^It is worth noting that asymptotically de Sitter solutions with squashed S'^ boundary metric, known as 
Taub-NUT de Sitter space, exist in Einstein gravity with a positive cosmological constant, see for example 
[45] . It would be interesting to explore the existence of such solutions in the bulk Vasiliev gravity. 
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3.1 Squashed geometry 

The metric of the squashed is given by an 5*^ fibered over an S'^ base space: 

ds'^ = J (^de^ + cos^ + ^^1^^ {dij + sin ^^0)^^ , (3.1) 

with V ~ + 47r, G [-f , f] and </) ~ </) + 27r. The squashed 5'=' is a homogenous yet 
anisotropic space. When a = the geometry becomes the round metric on expressed as 
a Hopf fibration. The Ricci scalar of the squashed is given by 

(3.2) 

r^(l + a) 

and becomes negative for a G (—1, —3/4). The volume is 

V=^=. (3.3) 
V 1 + « 

In what follows we set r = 1 unless otherwise specified. 

The eigenvalues of the conformal Laplacian on the squashed sphere are given by [46]: 



An,g 



+ a{n — 1 — 2g)' 



4(l + a)_ 



g = 0, l,...n-l ,n = 1,2, ... , (3.4) 



with multiplicity n. The values of a which have vanishing eigenvalues are all less than 
a = —3/4 and in fact accumulate rapidly as a — t- —1. Based on our calculation of the 
previous section we might expect wild oscillation for a G (—1, —3/4). 

The determinant of the conformal Laplacian for the squashed sphere has been studied 
using zeta function regularization in [47, 48], but it will be convenient for us to use a more 
numerical approach. 



3.2 A New Regulator 

We could attempt to study the squashed sphere determinant using a simple hard eigenmode 
cutoff as in equation (2.7) above, but this method is too clumsy to practically deal with the 
less symmetric geometry. We will instead use a heat-kernel type regulator [49], which for a 
general set of eigenvalues A, defines 

2 log Zj.r^ue = -2S,t -NJ2 y e-*^' . (3.5) 
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The justification for tliis definition is that for ^ 1 we have 

- j je''^" = -r(0, Xns) = log {Xnse"') + 0(A„£), (3.6) 



while for Ae ^ 1 we have 

" 1 



-r(0,Ae) = -e 



-Xe 



o 



Xe \{Xe) 



(3.7) 



Thus this proposal is equivalent to the determinant for modes whose energies are less than 
a "soft" cutoff and it cuts off the sum exponentially above the cutoff. The factor of 

e"^ is a wave function renormalization, similar to a in the previous section. 

It is convenient to split the integral over t into two pieces, which for the squashed sphere 
have the form 

2 log Z finite = detuv + det/R , (3.8) 



with 



n=0 q=0 



e 

oo n— 1 ,, oo n— 1 



det,^^-iVj]J]n / -e-*V, = _ivJ]J]nr(0,A„,,5) . (3.10) 

n=0 9=0 "^"^ n=0 (j=0 

Here 5 is some small but 0(1) number, say .01. The sum in det/^j is easy to perform 
numerically, one needs to sum up to about n = 1/5 to get good accuracy, detuv requires a 
little more care, since the e dependence needs to be treated analytically in order to extract 
the counterterms. Our method is to apply the Euler-Maclaurin formula^^ 

Y: m ~ / dnfin) + I (/(a) + /(6)) + {f^''-'\b) - f^''-'\a)) (3.11) 

n=a fc=l ^ '' 

to each sum in det[/y, before doing the integral over t. It is straightforward to see that the 
corrections to the formula are of increasing power in t, so if 5 is small the integral is a good 
approximation to the sum. In our calculations we included terms up to A; = 2 for both sums. 
We can expand the result in powers of t and then perform the integral; clearly only the lowest 
powers of t can produce UV divergences so we can extract the divergences analytically and 

^"'^In this formula B2k are the Bernouhi nunibers, and the series on the right-hand side is asymptotic. 
Truncating it at some finite k, the error is of order dn\f^'^'^\n)\. 
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Figure 2: A plot of l^'p with = 2 as a function of p, with p related to the squashing a as 
a = e^P- 1. 

subtract them with an appropriate choice of A and B. Conveniently there is no quadratic 
divergence, so the wave function renormahzation is apparently automatically handled by 
this regulator. Finally we can take the limit e — )■ and add this analytic expression to our 
numerical result for det//j to get a plot of \Zfinite\'^ = I'^hhI'^ over a wide range of values of 
a. We show the result in figure 2. 

There are some things to notice about this picture: 

• It is peaked at the round sphere a = 0. Apparently the universe prefers being round 
to being squashed, at least when we condition all of the other sources to be zero. The 
negative-definiteness of the stress tensor two-point function in fact shows that this 
peak is perturbatively stable under arbitrary metric deformations, with the only fiat 
directions corresponding to three-dimensional diffeomorphisms and Weyl rescalings.^^ 
The peak gets sharper as becomes large. 

• It is exponentially suppressed both for large a and for a — )■ — 1; there does not appear 
to be any instability for a < — |. If one zooms in on this region oscillation is indeed 
seen as in the mass deformation case, but the peaks are exponentially damped rather 
than exponentially growing. 

^■^The simplest way to see this negative definiteness is to recall that this partition function is the inverse of 
the 0{N) model partition function, which has a reflection-positive path integral and a stress-tensor two-point 
function which is positive definite. 
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• In fact we can argue that the round metric is a global maximum of the partition 
function over metrics with positive, but not necessarily constant, Ricci scalar. The 
reason is that the negative definiteness of the stress tensor two-point function is true 
for any metric with such a Ricci scalar, so the partition function is real and weakly 
concave. Any maximum must thus be a global maximum. ^'^ This result is consistent 
with [50] from the mathematics literature. Unfortunately we cannot decisively rule out 
other maxima with non-positive Ricci scalar which are higher than the round sphere. 
It would be very interesting to explore this possibility. 

4 The S'^ X partition function of the A^)' model 

We now turn to S"^ x S^}^ Because of the topology the singlet constraint is now nontrivial; 
we need to deal with the Chern-Simons coupling. The idea is that weakly gauging the Sp{N) 
symmetry will dynamically enforce the singlet constraint [54], and that in the limit A; — )■ oo 
it will have no other affect on the dynamics. We thus want to compute 

Zfinite X S^] = lim e-^^' [ VAVx e-^c:s[A]-5,„.e[x,A,9] _ (4^^) 

Share depends on the gauge field Ai through covariant derivatives, while the Chern-Simons 
action Scs is independent of the metric: 

Scs = ^^j d?x e^^'^Tr (^A,{d,Ak - d^^A^) + ^A,[A„ A^]^ . (4.2) 

In the limit A; — )■ oo we can evaluate the Chern-Simons path integral semiclassically by 
looking for flat connections. On S*^ the only such connections are gauge equivalent to Ai = 0, 
and the two sectors decouple. On 5*^ x 5*^ there are flat connections which have nontrivial 
holonomy around 5*^, and which have zero Chern-Simons action. These connections do not 
decouple in the A; — ?■ oo limit and must be integrated over explicitly. Our calculations in this 
section make extensive use of the tools developed in [55, 56, 57, 58, 59], but we will keep the 
discussion reasonably self-contained. 

^■^The condition on the Ricci scalar is necessary to ensure that in the 0{N) model the path integral can 
be done on a real integration contour, which is needed to conclude that its two-point function is positive. It 
really is necessary, since for a < — | there are indeed other maxima and the concavity breaks down - in this 
case however these other maxima are only local and not global. 

^^We acknowledge here soon to appear related work [51], where it is shown that in three-dimensional bulk 
de Sitter gravity there is a result similar to our main result of this section; the wave function with boundary 
geometry diverges at large conformal structure. 
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It is very convenient at this point to slightly modify the theory under discussion from 
an Sp{N) theory of Grassmans to a U{N) theory.^'' x° becomes a complex grassman field 
and instead of bilinears like x^^afeX* the action and currents now have the form x'^X- It is 
very easy to go back and do all of our previous calculations again in this new model, the 
only difference is that the power of the determinant in equations (2.3), (2.5) is instead of 
N/2. No qualitative results are altered. What this modification buys us is that we may now 
do a f/ (A^) gauge transformation to diagonalize the holonomy matrix for the Chern-Simons 
saddle point. Say we work on S"^ x with coordinates such that the metric is 



(4.3) 



We will take A to have periodicity one, so the parameter /3 measures the relative size of to 
S"^. We will sometimes refer to its inverse as the temperature T, although the interpretation 
here has nothing to do with thermality. In conventions where the covariant derivative is 
Di = di + Ai, we may choose a gauge where any flat connection has the form 



A, 



\ 



-lai 








-iCi2 








—iaN 



\ 



(4.4) 



This gauge fixing costs a Fadeev/Poppov determinant, and there is a remaining integral over 
the eigenvalues a^. Integrating out the Grassman fields one finds [54, 57]: 



■'finite 



1 



an) exp 



g sm 



ct 



ndet(-(5. + An)(5^ + ^n) + ^) • 



(4.5) 



^'^This theory is related to the bosonic U{N) model by iV — > —N so it is sometimes whimsically called the 
U{~N) model. Group theoretically we have that U{—N) — U{N) and 0{—N) — Sp{N) with symmetric 
and anti-symmetric representations switched [52, 53]. 
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To calculate the determinants one needs the eigenvalues^^ 

\n,^ = + + /^''("^ + 2^'^)' n e Z, £ = 0, 1, 2, . . . . (4.6) 

The determinants can be regularized and renormalized using a hard cutoff as in equation 
(2.7) above, the essential steps are described in section four of [57]. The result is 



^finite 



where 



^ / \{{da^)e^^ J^logsin^ - 2 i-;^5(e-^-) cos(ma 

" n<m ^ ' m=l n=l 



•7) 



/ \ i I *^ /a \ 

zs{x)=x-^ ^^_ . (4.8) 



This differs from the bosonic U{N) result of [54] only in the sign of the second term in the 
exponent. 

At finite it is difficult to learn more from this expression, but at large one can do 
the integral over the eigenvalues ctj using semiclassical techniques. This was done in the 
AdS case by [54], whose work we carry over almost verbatim. There are two interesting 
regimes. For temperatures that are 0{N^) the first term in the exponent is 0{N'^) while 
the second is only 0{N), so the saddle point is a small perturbation from what it would 
be with no matter present. Since the first term is repulsive, the eigenvalues will be close to 
uniformly distributed. Once the temperature becomes of order \/N however the second term 
also becomes 0{N'^). This is because for large T we have zs{e~^'^)) ~ The eigenvalue 
distribution then starts to be dominated by the second term, which is a potential that pushes 
the eigenvalues to tt. We now study this more quantitatively in both cases. 

4.1 Case 1: T <^ ^/N 

We begin by considering 'low' temperatures T <C V^. Since in this case the eigenvalues are 
roughly uniformly distributed it is convenient [55] to introduce an eigenvalue distribution 
function p{a) which counts the density of eigenvalues around some particular value a. We 



^^Here we have assumed periodic boundary conditions for the fermions. Had we chosen antiperiodic 
boundary conditions we would have gotten the same resuh except with all a.; shifted by tt. This shift drops 
out of equation (4.7) below since it can be absorbed into a change of variables in the integral over the a^. The 
wave function is therefore independent of which sign choice we make for the fermion boundary conditions. 
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then have: 



log Z finite = J dad(3p{a)p{l3) log 



sm 



a — (3 



r °° 1 

- 2N / dap{a) ^ 

^ m=l 



zsix ) coslma ) . 



(4.9) 

The density of eigenvalues p{a) in the absence of the U{N) vector matter fields is simply 
p{a) = This leads to the following A^^ contribution: 



\ogZj,r,ue[T] = -NHog2 . 



(4.10) 



One can compute the 1/A^ sub leading correction to the density, which we denote as p/N, 
such that: 

P(«) = + ^P(«) • (4.11) 
Following the discussion of [59] we find the saddle point equation to be satisfied by p{a): 

' a — P 



P / d(3p{(3) cot 



-2 ^ zsix"^) sin(ma) 



m=l 



which is solved by: 



oo ^ 



(4.12) 



(4.13) 



m=l 



Here P means principal value. Thus the partition function receives a correction given by: 



sm 



6 log Zfinite[T] = J dap{a) ^ / d/3p(/3) log 

cos(ria) cos(m/3) 



a — (3 



oo 



zs{x jcoslma) 



m=l 



/ dad(3- 



n 



n 



log 



sm 



a — (3 



m=l 



(4.14) 



m 



This is the only piece of Zfinite[T] that depends on temperature and thus the geometric data 
of S"^ X S^. The first piece is somewhat subtle to evaluate. First we note that it vanishes for 
m ^ n. After some work we find: 



dadf3 cos(ma) cos(n/3) log 



sm 



a — (3 



0. 

m 



m,n • 



(4.15) 
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Thus, putting everything together: 



(4.17) 



Due to a subtle cancelation in the signs, the above term in fact agrees with the answer for 
AdS4 [54] in the low temperature regime. It is the partition function of a free gas of higher 
spin particles. 

4.2 Case 2: T > ViV 

As mentioned above, the above analysis breaks down at sufficiently large temperatures since 
in that case the ~ A^^ and ~ A^ terms begin to compete. In the bosonic U{N) model one 
finds a Gross- Witten [59] type transition for which the density of eigenvalues strongly favors 
that they all vanish. Interestingly, due to the sign flip in front of A^ in (4.7) the potential at 
very large temperatures is V{a) ~ — cos a and forces the eigenvalues to be near a ~ tt. We 
can then approximate (4.7) by 



where . . . are terms that are either subleading in A^ or 1/T or are temperature independent. 
As in the low temperature case, the result is similar to the bosonic U (N) model studied in 
[54] in that the partition function grows with increasing T. In flgure 3 we plot 6 log Z finite[T] 
for small and large T. We again make a few comments about this result: 

• There does not seem to be even a local maximum for any flnite value of T. The prob- 
ability distribution apparently prefers arbitrarily small radii for 5*^. This divergence is 
more subtle to interpret than that of the scalar; since the topology is not something 
that can be embedded into de Sitter space there is no particular reason to view it as 
an instability. Rather it is a statement about a strange type of "big bang" where the 
universe is created with x S"^ topology. Perhaps surprisingly, we will show in section 
5 that it is actually possible to reproduce this type of divergence in standard Einstein 
gravity with positive Abuik- 
^''To obtain the result (4.17) we have used that: 



logZf,^ue[T] = 3mNT^ + ... , 



(4.18) 




(4.16) 



for non-zero integer m. 
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Figure 3: Left figure: Zfinite[T] for T < ViV. Right figure: Zfinite[T] for T > \/]V with iV = 1. 
For larger A'^ the curves become increasingly steep. 

• The divergence depends crucially on the topology - the squashed sphere also locally has 
the isometry of x S"^ but small S^, which corresponded to a — )■ —1, was exponentially 
suppressed. 

• In the bosonic U{N) model the result at large temperature is [54] log Z finite = 4^(3) A^T^. 
This is clearly not related to (4.18) by analytic continuation — t- — A^; the reason is 
that the eigenvalues cluster at vr instead of 0. We view this as a nice demonstration that 
dS and AdS are not in general related by analytic continuation at the nonperturbative 
level, even though they are in perturbation theory in many cases [5, 19, 23]. 

5 Comparison with Einstein Gravity 

In this section we review the Hartle-Hawking wave function for a free massive scalar in dS4 
with boundary geometry and for gravity with boundary geometry S"^ x S^. 

5.1 Free Scalar 

The wave function of a free massive scalar with future geometry can be lifted from the 
AdS calculation done in appendix C of [19], along with the analytic continuations X = 
and iads = iC-ds- This is in coordinates with metric 

ds'^ = -dr^ + cosh^ rdnl , (5.1) 
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where we have set ids = 1- Speciahzing the general result of [19] to d = 3 and A = 2, and 
also assuming that only the zero mode is turned on, we have 

log^'[r,0] = ivr^cosh^r^o [-1 + ^(sinhr)"^ + . . .] . (5.2) 

To compare to a CFT calculation we define ao = 0coshr, after which we see that 

= e""'^' . (5.3) 

Changing to the a basis we have 

= . (5.4) 

The points to notice are that the divergent counterterm is pure phase and the wave function 
is square-normalizeable. There is no sign of the divergence we found in the Sp{N) model as 
a function of the mass deformation. 



5.2 Einstein Gravity with Positive Cosmological Constant and "Thermal" Bound- 
ary Conditions 

In this section we present a Ahui^ > version of the semiclassical Einstein gravity calculation 
[60, 61] of the wave function with S'^^^ x S'^ boundary conditions We will specialize to d = 3 
at the end. We are interested in compact complex solutions of the Lorentzian action 



Sr 



IGnG 



j d^+^xy/^{R-d{d-l)) + 2 J d'^Xy/^K 



(5.5) 



which have a single boundary with topology S x S and induced boundary metric 



jijdx'^dx^ = r1 



pde 



(5.6) 



We have chosen the cosmological constant here so that the dS radius is one. The action 
evaluated on such a solution is 



Sr 



1 



SttG 



d I d'^^^x^ + 2 / d^x^K 



(5.7) 
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The only solutions with SO{d) x SO {2) symmetry are analytic continuations of dS/Schwarzchild: 

dp 



ds' 



fir] 



+ f{r)dX' + r'dn 



d-2 5 



with 



/(r) = — 1 + ar 



2-d 



(5.8) 



(5.9) 



The coordinate r runs over a range r G (ro, Tc), where for there to be no boundary at tq we 
must choose ro to be a root of f{r). The coordinate A we take to be periodic with period 
Aq. The trace of the extrinsic curvature at the boundary is 



K 



(5.10) 



and the action evaluates to 



iSr 



SttG 



d 



d-2 



id - lyr^rc' - 1) - ^rt' + 2 



(5.11) 



To match the periodicity of the at the boundary, we must have 



An 



1 



1 



(5.12) 



where . . . indicates terms that fall off faster as Tc — )■ oo. Inserting this into the action, we 
find 



iSr 



i^d-i(^ 1 



SttG 



1 



{d - lyt^ri - 1) + -rt\l + rt) + 0(l/r,) 



(5.13) 



Note that all divergent terms as Tc — ?■ oo are pure phase, as argued for generally above 
equation (2.5). In odd dimensions the square root in the denominator just shifts around the 
coefficients of the divergences, while in even dimensions it also contributes a phase to the 
finite piece. 

Finally we observe that ro and Ao are related by the requirement that the (complex) 
geometry be regular at ro. There are two ways for this to happen. The first is if ro=0; in 
this case a = and the geometry closes off because the S'^~^ shrinks to zero size. We may 
then set Aq at will to match /3 as in equation (5.12). This geometry is a quotient of pure 
de Sitter space, and its wave function ^ = e*'^^ is pure phase. The other option is for the 
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to shrink to zero size at finite radius of the S'^ ^, but in a way that the geometry caps 
smoothly. This requires 



' 

or equivalently 



Note that there are four possible complex solutions. Intuitively the relative sign arises 
from whether tq is placed at the black hole or cosmological horizon of the continued de 
Sitter/Schwarzchild, while the overall sign is a choice of branch in how we continue the regu- 
larity condition away from real Euclidean geometries. Without some sort of non-perturbative 
information about the Einstein gravity path integral we have no way of deciding how many 
of these solutions should contribute to the semiclassical limit of the wave function. The 
most conservative thing to do is to acknowledge that any of them give a valid solution of the 
Wheeler deWitt equation, and study each in turn. 

We are especially interested in the limits of high and low temperature. For low temper- 
ature we have tq — ?■ which interestingly is the Nariai value of tq. The finite part of 
the action at low temperature is then 

^%..e = ^^(±) — (— • (5.16) 



This is a pure phase in the wave function that oscillates as /3 — oo. The behavior at 

d-2] 

4lT 



high temperature is more interesting, there are two types of solution: tq — > i^i^L^ or 



To — > ±-^. In the first type the action vanishes as /3 — > 0. In the second type the action 
diverg GS, clS db phase in the wave function for even dimension and as a diverging real part in 
odd dimension. ^'^ In particular for d = 3, the behavior is 

167r^ 

finite = ^27GP^ ' (5-17) 

For the upper choice of sign this is exactly the scaling with G and /3 that we found in the 
fermionic U{N) model in equation (4.18)! Apparently the divergence of the wave function 
at large temperature in the dual field theory is capturing something about de Sitter space 
that is present even in Einstein gravity. It is interesting to note that the saddle points which 

^*^The divergence found by Castro and Maloney [51] for d = 2 does not contradict this result because the 
solutions we are considering are only valid for d > 3. 
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control this limit are dS/Schwarzchild black holes with large and purely imaginary mass. 
6 More on the Scalar Divergence 

In this section we return to the divergence at large negative ctq we found in section 2.1. We 
first discuss the extent to which it can be interpreted as a non-normalizeability of the wave 
function. We then study the partition functions of the critical 0{N) and Sp{N) theories 
on as a function of a constant mass deformation. We'll see that in the Sp{N) case it is 
actually impossible to define the critical theory on nonperturbatively in A^. Finally we 
will return to the question of the bulk interpretation of the scalar sector, concluding that the 
most reasonable interpretation is that it indicates a non-perturbative instability of Vasiliev's 
theory in de Sitter space. 

6.1 Normalizeability? 

Because we are computing conditional probabilities where the sources for higher-spin fields 
are set to zero, it is possible that the exponential divergence we found in section 2.1 does 
not actually indicate non-normalizeability of the wave function. For example the function 
, is integrable near x = y = 0, but if we set y = and try to perform the x integral 
we will find a spurious divergence. We can fix this by setting x to be non-zero and then 
integrating over y in some finite region around y = 0. The resulting distribution for x will 
then be integrable. In this subsection we will suggest an argument that integrating over 
any finite number of the zero modes of the higher spin sources will not resolve the apparent 
non-normalizeability of the wave function in the scalar zero mode, but to make this argument 
precise the higher-spin sources need to be studied in much more detail than we provide here. 

In the presence of a constant mass deformation the higher-spin currents can be modified 
in such a way that they are still covariantly conserved. For example the stress tensor (1.3) 
picks up an extra term ^crx"x''- More generally the spin s current has the heuristic form 

J^'\,...,. = ^at {x'd,, . . . d.a" + ...+ g,,,, . . . g,,_,,^a^/W) , (6.1) 

The indices are symmetrized and have traces subtracted in the appropriate way. We can 
turn on linearized sources for these operators in the action 

SScFT = J2j d'x^9h^^- '^f'\^,„^^ . (6.2) 
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To Gaussian order the source dependent piece of the partition function is 

At finite a the two-point function of J'-*-* in this expression is calculated using massive 
propagators. Various quadratic counterterms in h^^_^^^^ are needed to to make the correlator 
finite which we will not be explicit about. Evaluating this integral on 5''^ is difficult, it 
requires higher-spin spherical harmonics. On however it is easy to see that the integral 
is proportional to a power of a. Performing the Gaussian integral over h^^ will bring 
down the square root of the inverse of this power of a into the integration measure for a, 
but it will never be able to beat an exponential divergence of the type found in section 2.1. 
Unfortunately although the answer on M? is a good approximation to the calculation 
in most of the a plane at large cr, it fails precisely on the negative real a axis where the 
divergence happens. The full calculation is beyond the scope of this work, but we note 
that generically one encounters powers of a times sums of the form 

oo ^ 

This sum can be performed analytically, and its inverse does not approach zero at large 
negative a. Instead it oscillates with constant amplitude. To the extent that this sum is a 
good model for the full set of terms that appear we see that integrating over /i/^i.../^^ will not 
affect the non-normalizeability of the integral over a. 

6.2 The Partition Function of the critical 0{N) Theory on the Three Sphere 

As explained in appendix A, one expects the wave function in the field basis to be computed 
by the interacting critical Sp{N) theory reached in the IR of a "double trace" flow from the 
free Sp{N) theory. This theory ends up being subtly inconsistent, so as a warmup in this 
section we compute the partition function on as a function of a constant mass deformation 
of the bosonic critical 0{N) model conjectured to be dual to higher spin gravity in AdS4. 
The tools we develop for this more familiar case we will then apply to the Sp{N) theory in 
the following subsection. 
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We first define 



^ free 



(T = e 



^NI[a] 



l?0exp 



and 



■'crit 



Vcpexp 



(6.5) 



(6.6) 



The 0's have vector indices running from 1 to that are contracted in the obvious way and 
/ G C is some constant. The function /(a) is defined as the right hand side of equation 
(2.13) divided by A^. Note the free 0{N) partition function is the inverse of the free Sp{N) 
partition function. As discussed in more detail in the appendix the critical partition function 
is a wave function (here a radial one since we are in AdS) in the field basis. These partition 
functions are related by 



Zcrit [c] 



'det 



-A^ 



27r/ 



e 32 



Va exp 



A^ / rffi. 



V2/ 4 



era 



27r2 



o- 



, (6.7) 



which is a special case of (A. 8) in the appendix. At large A^ this integral can be evaluated 
semiclassically by looking for a saddle point ac of the integral over a. When a is constant, 
the saddle point equation is be algebraic (see for example the appendix of [44]): 



47rcro = a/1 - 4crc cot (^^Vl - 4(7, 



(6.^ 



For real ao this equation has infinitely many solutions, one greater than —3/4 and the rest less 
than it. The one which is greater than —3/4 approaches —3/4 asao — oo and approaches 
+00 as 5^0 ~^ +00. To determine which is the relevant one to use in the semiclassical 
approximation we need to study the integration contour for cr. For / real and positive we 
must take the a contour to come in from a wedge of angle tt/2 centered on the negative 
imaginary axis and go out through a wedge of angle 7r/2 centered on the positive imaginary 
axis. For example we can just take it to run along the imaginary axis. As / becomes large 
however the behavior of the integrand is controlled by J(cr) for a very long time before the 

2 

term |y takes over. For the integral to have a limit as / — )■ cxd it is necessary that we 
can deform the original integration contour smoothly into one that remains convergent as 
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In this subsection only is a boundary field, not a bulk field. 
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Figure 4: Important features of the a plane, and a plot of our large result for 



exp 



as a function of a. 



^logZeritp]-^ 

/ —7- oo.^" The asymptotic behavior for /(cr) away from the negative real axis is 

I{(J) ~ -^fxi , (6.9) 

D 

SO in order for the the limit / — )■ oo to exist the defining contour must be deformable to 
one that enters from the region — tt < 9 < — 7r/3 of the sigma plane and exits through the 
region tx /?> < 9 < ti. We shade these regions in blue in figure 4, which illustrates the relevant 
features of the a plane. The singularities of the integrand, which are the values of a where 
the determinant vanishes, are marked as x's. We are not allowed to deform the integration 
contour across them. The saddle point with > —3/4 is marked in red for a typical value 
of aQ. The lesson however is that it is clear from the figure that the defining contour can 
always be deformed into the red stationary phase contour that passes through the red saddle 
point, so we can approximate Zcrit at large N as the contribution of this saddle point only. 
We plot our result for the finite part of Zcrit in figure 4. For generic Oq the saddle point 
equation must be solved numerically, but in asymptotic limits it is not hard to see that 

^ log Zcruia] - - ^ (log4 - ^) - vr^?^ ? (6.11) 

logZcritlcr] - '—a^ ^ —a - - log — a ^ -oo . (6.12) 



N ^ '16 8 2 ^ V 7r2? 
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The partition function will still be divergent in this limit because of the local term in front of the integral 
in equation (6.7), but we want to make sure that all divergence has been isolated in this term. 
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Finally we observe that when a = 0, the relevant saddle point is cTc = 0. This is consistent 
with the perturbation theory usually used to study correlation functions in the critical 0{N) 
model; recall that in this perturbation theory at large N one sums over "cactus" diagrams. 
For example the two-point function of 0^ at large is dominated by a set of diagrams shown 
in figure 5. Here the dashed line is a a propagator — y. These diagrams are a geometric 
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Figure 5: The first few diagrams that contribute at order N to the cf^ two point function. 

series which can be summed to see that 0^ has dimension two when the momentum transfer 
is low compared to /. Using these diagrammatics we can confirm that the one point function 
of a is indeed zero as long as we normal order 0^ in the action. 

6.3 The Partition Function of the Critical Sp{N) Theory on the Three Sphere 

We now attempt to study the critical Sp{N) theory using the same tools just applied to 
the 0{N) model. As we explain in detail in appendix A, we would like to interpret this 
partition function as the field basis representation of the Hartle-Hawking wave function. 
From equation (A. 8) we are interested in the integral 



ZrriilO'l 



'det 



2nf 



N / dn-, 



(6.13) 



where we have used the results of section A. 3. Here e is a real and positive length cutoff 
defined so that / = foe~^. We will get to the critical theory by taking e — )■ with /o fixed. 
The phase of /o is a free choice, we will consider various possibilities. The saddle point 
equation is 

Aniao = ^1 - 4ac cot (^^^/l^^Ia~^ , (6.14) 

so there is still a saddle point Uc which is zero when a = 0. As in the 0{N) model the 
"cactus" perturbation theory is based on this saddle point. There are again however lots of 
other saddle points; essentially there is one associated with each of the zeros of the functional 
determinant of the conformal Laplacian. In figure 6 we show the important features of the 
a plane analogous to those in the 0{N) model. 

There is a clear difference in that the allowed region in which the integration contour 
may approach infinity, again shaded blue, is much smaller. Unlike the AdS case, there are 
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Figure 6: Important features of the a plane for the dS case. 

no singularities in the blue region. This means that the integral on any contour that starts 
and ends at infinity in the blue region must be zero! The only way to get a nontrivial result 
on a contour that runs from infinity to infinity is for the contour to come from and/or go 
to infinity along the negative real axis, where the approximation (6.9) breaks down. For 
example if /o is real and negative then we can take the defining contour to run along the 

2 

real axis. The term |j in the exponent of the integrand will render the integral convergent, 
but the integral will be dominated by the region of very negative a. The limit / — )■ oo will 
not exist. It is not difficult to see that there is no choice of phase for /o where the defining 
contour can be deformed in such a way to produce a good limit as / — )■ oo. In particular there 
is no defining contour from infinity to infinity on which the integral can be semiclassically 
approximated by the cactus perturbation theory; the stationary phase contour which passes 
through the "cactus" saddle point runs from —3/4 to infinity and is indicated in figure 6. 
Thus we see that non-perturbatively in N the critical Sp{N) theory on 5*^ does not exist. 

6.4 Bulk Interpretation 

What are we to make of the scalar divergence and the non-existence of the critical Sp{N) 
theory at finite source? That both of these theories should have well-defined partition func- 
tions with finite sources turned on is equivalent to the statement that in the bulk there are 

^^One might be tempted to take the red contour in figure G as the defining contour for the critical Sp{N) 
theory. This would be giving up on the double trace fiow interpretation, and it would be unpleasantly 
ambiguous. Why should the path integral stop at cr = —3/4 and not some other value? We could try to 
make it sound less ambiguous by defining the end point as the least negative eigenvalue of the conformal 
Laplacian, but the path integral would then depend non-locally on the metric and calling it a field theory 
would be questionable. This type of metric dependence does not seem consistent with what we would expect 
from the bulk. We do not find options of this kind appealing. 
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limiting probability distributions for a and a which are well defined at late times. This 
is essentially the assumption that the late time behavior is asymptotically de Sitter in the 
Fefferman-Graham sense. We have uncovered a contradiction in this assumption for the 
higher-spin theory dual to the free Sp{N) model. The only way out is to conclude that this 
model does not in fact describe a theory which is asymptotically de Sitter at late times. 
What instead must happen is some sort of instability that allows a to find its way out to 
large negative values. Recall that in quantum mechanics non-normalizeability of the ground 
state wave function is the key indicator for an unstable system. For example a particle in an 
exponential potential has a ground state wave function which grows linearly with distance 
on the decreasing side of the exponential. Since the Hartle-Hawking wave function is in some 
sense supposed to be describing the equilibrium configuration of the theory, we conclude that 
in this equilibrium the field a tends to be far from its de Sitter invariant value. Of course 
since we were not able to decisively prove non-normalizeability for the technical reasons de- 
scribed in section 6.1, we can't completely rule out the possibility that once all of the sources 
for higher spin fields are integrated over there is a miraculous pressure to positive a that 
restores the de Sitter invariance. We have found no indication so far that this happens, and 
we find it simpler to envision that higher spin gravity in de Sitter space is unstable. 

7 Three- manifolds and the 'big bang'? 

Pure de Sitter space is known to be a solution of Vasiliev's higher spin gravity with only 
the metric turned on and all other fields turned off. It corresponds to the CFT living on 
some three-geometry with all sources turned off. Which three-geometry depends on how 
the boundary is approached. Some common choices are the standard metrics on 5*^, H^, 
M.^ and S*^ x M^. They correspond to the global, hyperbolic (relevant to the nucleation of a 
Coleman-de Luccia bubble), planar and static patch slicings of pure de Sitter space. 

In addition, one may consider quotients of the above geometries which are still solutions 
in Vasiliev gravity. Such quotients can make the de Sitter slicing compact. The simplest 
example is the tq = solution we studied in section 5, but there are many other options. 
For instance, we can quotient to a T^. This geometry shrinks to zero size in the infinite 
past, touching a point at X^. Another example corresponds to quotienting 'H'^ to a compact 
manifold. In this case the universe shrinks to a point at a finite proper time in the past. 
This is clear from the metric of de Sitter space with slices: 
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where dTil = [dp^ + cosh^p dVt^ is the canonical unit metric on The 'big bang' would 
occur at r = 0. 

From the perspective of the dual CFT, we are considering the Sp{N) theory on a compact 
quotient of Ti^. Given that we must impose a singlet constraint, it is rather complicated to 
actually compute the partition function as a function of the geometric data of the quotient. 
On the other hand, it is well known mathematically that such compact quotients have a rich 
topological structure which the integral over the space of flat connections will be sensitive 
to. It would be interesting to understand more deeply how the presence of a 'big bang' is 
connected to this topological structure. 

We should mention a criticism to the above. Vasiliev theories are non-local for distances 
smaller than the de Sitter length. This is due to an infinite sequence of higher derivative 
terms in the equations of motion and the fact that the de Sitter length is the only scale in the 
theory. The big bang singularity we are discussing is at scales far smaller than the de Sitter 
scale, so without the aid of a good gauge invariant observable we cannot conclude whether 
there indeed exists a pathology. On a more perturbative level however, these solutions have 
a large metric turned on and all other fields turned off, thus a probe like observer would 
indeed perceive a shrinking universe. Some insight may be gained from a linear analysis 
about one of these 'big bang' cosmologies. 

8 Outlook 

The nature of our work has been explorative. There are many directions in which we must 
proceed to flesh out our understanding of de Sitter holography, even in the context of the 
simple example we are studying. It would be interesting to understand how much of our 
picture is qualitatively similar to the case of more ordinary theories of Einstein gravity. An 
example of this is the divergence at small 5*^ relative to 5*^ discussed in section 4. 

An interesting question that naturally arises is that of topology. How are we to compare 
the wave function evaluated on different topologies? In situations with more complicated 
topologies the Chern-Simons piece of the action will play a significant role.'^^ A simple 
example that seems calculable is the partition function on a Lens space (i.e. a quotient of 
by the cyclic group of order p, Zp, which has finite fundamental group). Such Lens spaces 

^^An obvious guess for how to compare different topologies is to normalize the Chern-Simons integration 
measure on some manifold M as in [62]. In this normalization the pure Chern-Simons partition function on 

= S*^ X 5^ is Zcs [S'^ X S^] = 1 for any gauge group and k, whereas it vanishes at large k on an 
More generally this normalization would suggest that more sophisticated topologies are infinitely preferred 
as k ^ oo. We are unsure however if this is really the correct way to compare topologies. 
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provide a discretum of topologies that we can study the wave function on. More generally, 
inspired by [42], we might also envision studying the wave function on more general higher 
genus surfaces. 

We must also confront the issue of turning on sources for the higher spin currents. One 
approach would be to sort out the higher spin spherical harmonics needed to complete the 
argument of section 6.1 and decisively confirm the non-normalizeability of the wave function. 
Alternatively, given that the higher spin currents are all quadratic in the fields x", one expects 
schematically that the action for general sources is given by the bilocal expression [63, 64]: 



The challenge of course is to understand what we mean be the space of all B{x,y)^s and 



function whose coordinates are the sources of each current, but this formalism may simplify 
the treatment of the higher spin sources. 

Finally if the scalar divergence indeed indicates an instability of Vasiliev's theory in de 
Sitter space, it would be very interesting to understand the bulk mechanism by which this 
decay proceeds. In particular we would like to know what is the fate of Vasiliev's universe; 
where does it decay to? It has been recently argued [65] that the AdS version of Vasiliev's 
higher-spin theories can be realized as various limits of string theory; perhaps these decays 
connect higher-spin theories dynamically to the rest of the string landscape [71, 72, 73]. 
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how to relate the partition function Zcft 



J Vxexp{-Sbare[B,gij,x]) back to a wave 
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do not necessarily reflect the views of the John Templeton Foundation. 
A Double Trace Flows in dS/CFT 

In AdS/CFT it has long been understood [66] that there are two "natural" boundary condi- 
tions for a bulk scalar field, which are usually referred to as the "standard" and "alternate" 
quantizations. The partition functions for the two boundary conditions are related by a 
double-trace RG flow [67]. In dS/CFT the analogue of the standard quantization of a bulk 
scalar is the Hartle-Hawking wave function in field space, while the analogue of the alternate 
quantization is the same wave function projected onto a different basis of the Hilbert space 
of field configurations. In this appendix we develop this connection on the field theory side, 
showing in more detail how it arises from a double-trace RG flow. We will give arguments 
valid for a general large- CFT in d dimensions, which we will specialize to the Vasiliev case 
theory at the end. 

In AdS / CFT the observation that boundary conditions in the bulk are a statement about 
states in the radial bulk Hilbert space goes back at least to [68]; our discussion here builds 
primarily on [67] and [69]. These issues were also recently studied in the 0{N) model in 
[70], whose results we use to resolve an important technical subtlety. 

A.l Field Theory Double- Trace Flow 

Consider a CFT on M'^ that is deformed by a scalar operator O having dimension A < d/2: 

Zait[(j] = j i5Me-^^^^+^'^'""(")^(") . (A.l) 

For technical reasons we will also assume that A > max (|, ^^).'^'^ Say that the lowest 
dimension non-identity operator appearing in the OO operator product expansion is a scalar, 
which we denote (9^, and that there is some parameter A^ such that its scaling dimension 
is A2 = 2 A + 0(1/A^). We will normalize O so that its two-point function is 0{N)^^^ It is 
widely believed [68, 67] that at finite but large enough A^ the deformed theory 

Zfio^i^] = j PA^e^^^^^-^^'^'"^' (A.2) 

^^This inequality is actually saturated for the Sp{N) model; we comment on this below. 
^''^We have defined TV in a way that suggests a vector theory, but this section is also valid for matrix 
theories provided one replaces N — >■ N'^. 
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flows to an IR fixed point at which the dimension of O becomes A = d — A + 0{1/N). The 
scale below which the theory is near the IR fixed point is determined by the dimensionful 
parameter /, so by choosing / to be of order the cutoff length to the appropriate power we 
can directly study the IR fixed point. To be definite we will set 



Ao-d 



(A.3) 



with /o some dimensionless constant that is 0{N^). We can then study a scalar deformation 
of this fixed point by computing^^ 



^.tan[?] = e'^-[^l / VMexp 



-ScFT + pfe-'aO-^ I d'^xO^ 



(A.4) 



Here 7 = (i + A — A — A2 = 0{1/N), which ensures that a has dimension rf — A as expected 
from the dS/CFT dictionary, p is a dimensionless constant we introduce so that we may 
take a to be related to a canonically normalized bulk field in the standard way: 



Ne' 



i-A 



a 



(A.5) 



The value of p will of course depend on the CFT normalization of O. The counterterm action 
Set is included to ensure that we get the full bulk wave function including local terms. For 
the range of A we are studying there is only one counterterm that does not vanish as e — )■ 0, 
so we simply have 



5-, 



ct 



'-aNe' 



d-2A 



d'^xa^' 



(A.6) 



with a a tuneable dimensionless parameter we determine below by comparison with the bulk. 

Following [67] we observe that that Zstan and Zait are related by Hubbard-Stratonovich 
transformations. To simplify formulas we will from here on neglect 0(1/A^) corrections in 
powers of the cutoff; these can be restored at any point by dimensional analysis. We may 



^^In this equation and the previous one we should include a "bare source" proportional to O in the action, 
which we tune to ensure that the one point function of O is zero at the IR fixed point. The coefficient of 
this term is zero to leading order in A^, so rather than carry it around explicitly we instead declare that in 
what follows one should understand our a to differ from the "true" a by an additive constant that vanishes 
at large N. In equation (A. 14) below this translates into a small constant additive shift of a; this has no 
meaningful effect on the physics of the main text so we ignore it from here out. 

^^We will also assume that we can freely multiply Cs together using O x O = e'^^^^'^O^ . We should in 
principle also include the UV-divergent identity contribution to OO, but this gives only scheme-dependent 
factors which are independent of a and a and thus can be absorbed into the normalization of Zgtan and Zait 
in a trivial way. 
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then write: 



Z.y^-Pl+^Z'^''^^^ j Vae''^''<^-''''~^)z^u[a] • (A.8) 



ZstaJo-] = i/det 



The integration contours depend on the phase of / and are chosen to ensure convergence. 
It was pointed out in [67] that in the hmit e — > (which here just means / — )■ oo) the 
transformation (A.8) becomes a functional Legendre transformation up to a local countert- 
erm; the inverse transformation (A.7) becomes singular in a nontrivial way. In the following 
subsection we will interpret these transformations as a quantum mechanical change of basis. 

A. 2 Bulk Boundary Conditions 

In this section we recall some basic properties of the wave function of a massive scalar in dS 
space. For simplicity we will use the formalism of a scalar field in a fixed dS background, 
but we could also make the argument in Wheeler-deWitt language. We will use coordinates 
where 

(A.9) 



ds^ = — {-dT^ + dx") 



For a scalar action of the form 



bulk 



-d^<i)d^<j) + V{<j)) 



(A.IO) 



with T < 0, we can discretize the bulk path integral expression of the Bunch-Davies wave 
function to find 



vi>[0,r + 5] 



-T) 



exp 



d'^x 



V(t) 



- (5(-T)-(^+i) ( v(0) + y 



vi>[0,r] 

(A.ll) 



to linear order in 5. Expanding both sides in 5 one can derive the functional Schrodinger 
equation for We can easily see that canonical momentum conjugate to is 



n = -i- 



(A.12) 
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Our introduction of Set and p in the previous section allow a simple dS/CFT bulk-boundary 
dictionary: 

Zstanla] = ^VNae''-^, -e] . (A.13) 

Using this dictionary we can apply the transformation (A. 7) to the (A. 11) to find a bulk 
interpretation for Zait- The integral over a is Gaussian and can be done exactly - it imposes 
boundary conditions ^' 

- TdT(j) + i(pVo - a)0 = ip^a^-Tf-^ . (A. 14) 

Thus Zait is the same bulk path integral but with different boundary conditions. For real 
field configurations the left hand side of (A. 14) should be real. As we noted before the value 
of p will depend on the choice of normalization of (9 - we here set the phase of {00)ait such 
that this reality requires cr to be real. We then see that p and (p^/o — a) must be purely 
imaginary. 

With this phase choice we can then interpret Zait as the Hartle-Hawking wave function 
projected onto eigenstates of the hermitian operator 

S = ' (-(-^)'n + ^(pVo - «)0) . (A.15) 

pV A 

Note that if we had not imposed a lower bound on A, there could have been additional 
terms in Set and the boundary conditions could be more complicated. 

^^As mentioned in a previous footnote, because the Sp{N) model saturates the allowed region for A it is 
in principle possible to have an additional term in Set proportional to a^. This would lead to an additional 
term proportional to N'^cf)'^ on the left hand side of (A. 14). This term would be of importance to us since 
we are interested in field configurations where (j) y/N. The non-existence of such a term can be inferred 
from the work of [70], who showed order by order in that the critical and free 0{N) models are related 
in the same way as the bulk correlation functions built from propagators with Dirichlet boundary conditions 
or (A. 14) respectively.* This agreement would be spoiled by a quadratic term in the boundary conditions. 
The argument of [70] used special properties of the 0{N) model, which carry over straightforwardly to the 
Sp{N) model. It is worth mentioning that the interpretation of this term, according to the general proposal 
for multitrace boundary conditions in [G8], is that were it included it would deform the free Sp{N) model 
by the approximately marginal operator (x'^^abX^)^ ■ 

*Actually [70] don't explicitly mention (A. 14), the detailed equivalence of the "alternate" propagators 
they use to (A. 14) for the AdS case was shown earlier in [69]. 
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A. 3 Explicit Parameters at large N 

At large N we can determine p and a by comparison of {00)stan to a free bulk scalar 
computation of \t'. The needed bulk result is [19] 



^[0,T] = exp 



(A.16) 



where c = ^-^^ and . . . means terms that are higher order in — T. A is related 

d-A r(A-d/2) 

to the bulk mass in the usual dS/CFT way: 

A = ^ + - 4m2 . (A.17) 

We can compare this with the standard large N result [67] 



1 



/o /(^ 



(O.0.').tan = + ( - - i^ike)'-'^ + ...); (A.18) 

if we define 



|2A 

then we have [67] 



{0{x)0{y)U = (A.19) 
\x — y\ 



By matching the second derivative of Zstan to (A.16) we then find 

a-p'fo = tA + 0{l/N) , (A.21) 

and 

Finally for use in the main text we note that in the Sp{N) model we have = 3, A = 1, 
and if we define 

O = J'^'^ = -^-Qatx'^x' , (A.23) 
then C = -3^. This gives = _i/i6 + 0(1/A) and a - fop^ = i + 0(1/A), so the 
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boundary conditions become^^ 

TdT<p + = ^ViVTV . (A. 24) 

In fact for the Sp{N) model the 0{1/N) corrections to p and a actually are zero; this follows 
from the argument of [70] , and uses special properties of the Sp{N) model. This cancellation 
is not necessary for the general interpretation of the alternate quantization, but it makes 
things more convenient in the main text so we will use it. 
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